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ABSTRACT 

We present the results of an intensive ground-based spectrophotometric monitoring 
campaign of the Seyfert galaxy NGC 4151 for a period of over two months, with a typical 
temporal resolution of one day. Light curves for four optical continuum bands and the Ha 
and H/3 emission lines are given. During the monitoring period, the continuum at 6925 A 
varied by ~17% while the continuum at 4600 A varied by ~35%, with larger variations 
in the near UV. The wavelength dependence of the variation amplitude also extends into 
the far UV. The dependence in the 2700-7200 A range can be explained by the different 
relative starlight contributions at different wavelengths, but the large variability at 1275 A 
cannot be explained in this way. The continuum variability timescale is of order 13 days 
and is similar at all optical wavelength bands. No evidence for a time lag between the 
optical continuum and the UV continuum and emission lines was found. The Ha emission 
line flux varied by about 12% with a gradual rise throughout the campaign. Its cross 
correlation with the continuum light curve gives a lag of — 2 days. The variations in 
the H/3 emission line flux are about 30% and lag the continuum by 0—3 days. This is in 
contrast to past results where a time lag of 9±2 days was found for both emission lines. 
This may be due to a different variability timescale of the ionizing continuum, or to a real 
change in the BLR gas distribution in the 5.5 years interval between the two campaigns. 

Subject headings: galaxies: individual (NGC 4151) — galaxies: active — galaxies: Seyfert 
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1. Introduction 

The Seyfert 1 galaxy NGC 4151 is one of the best- 
studied active galactic nuclei (AGN) due to its bright- 
ness and variability properties. It has been studied 
at many wavelengths and its characteristics are well 
known (e.g., Peterson 1988). Several monitoring cam- 
paigns have shown variability timescales from a few 
hours in the hard X-rays (Yaqoob et al. 1993), to a 
few days in the ultraviolet (e.g., Clavel et al. 1990) 
and the optical (e.g., Maoz et al. 1991), and several 
months in the IR (Prestwich, Wright, & Joseph 1992). 

NGC 4151 was selected by the AGN Watch con- 
sortium as a prime target for an intensive spectro- 
scopic multiwavelength monitoring campaign. The 
campaign took place for two weeks in 1993 Decem- 
ber using the CGRO, ASCA, ROSAT, and lUE 
satellites, and many ground-based telescopes. The 
Crenshaw et al. paper in this issue (hereafter Pa- 
per I) describes the UV results from lUE. This 
paper describes the observations and results of the 
ground-based optical campaign. The high-energy re- 
sults {GEO, ASCA, ROSAT) are described by War- 
wick et al. in this issue (Paper III), and a multiwave- 
length comparison is given by Edelson et al. (Paper 
IV). 

A primary goal of AGN monitoring has been to de- 
termine the size of the broad-line region (BLR; see 
Peterson 1993 for a review). Of the many ground- 
based variability studies of NGC 4151, the most in- 
tensive have been those of Antonucci & Cohen (1983) 
and Maoz etal. (1991). Antonucci & Cohen monitored 
NGC 4151 at approximately monthly intervals for over 
a year. They found that the continuum and broad lines 
varied on timescales shorter than their temporal res- 
olution, and deduced a BLR radius of less than ~30 
It-days. For these data Gaskell & Sparke (1986) re- 
ported an H/3 lag of — 7 days and H7 lag of 5 — 9 
days. These observations were also analyzed by Peter- 
son & Cota (1988) who combined them with their own 
data and arrived at a BLR size of 6±4 It-days. The 
Maoz et al. (1991) monitoring campaign lasted for a 
period of 8 months with a mean sampling interval of 
4 days. Cross correlation, deconvolution, and model- 
ing applied to the data indicated a BLR size of 9±2 
It-days. The lUE monitoring campaign described by 
Clavel et al. (1990) yielded a characteristic timescale 
of 4 ± 3 days from the response of the C iv A1549 and 
Mgii A2798 emission lines to the UV continuum. 

In this paper we present results of two months of 
monitoring of NGC 4151 with a time resolution of 
about one day. In §2 we describe the ground-based 
observations together with their reduction and calibra- 
tion, and present light curves for the continuum and 
the Ha and H/? emission lines. In §3 we carry out a 
time-series analysis of the data and briefly discuss the 
results. A summary is given in §4. 

2. Observations 

2.1. Data and Reduction 

The main effort of the ground-based monitoring 
campaign was carried out at the Wise Observatory in 



Israel, using the Im telescope. The last two hours of 
each night, for a period of over two months starting on 
1993 November 14, were dedicated to the monitoring 
of NGC 4151. Spectroscopic observations were per- 
formed with the Faint Object Spectrographic Camera 
(Kaspi et al. 1995), using a 10"-wide slit and a 600 
line/mm grism, giving a dispersion of 2 A per pixel. 
The spectral resolution was determined by the seeing, 
which was 2 to 3 arcsec, and combined with the in- 
strument spatial scale of 0.9 arcsec/pixel the spectral 
resolution was ^5 A. Two different slits were used ev- 
ery night, each located at a different position in the 
telescope's focal plane. One slit produced spectra in 
the range ~ 4150-6050 A (hereafter the "B-side") and 
the other covered the range ~ 5060-6990 A (the "R- 
side"). With this setup we were able to monitor both 
the Ha and the H/3 emission lines, while a large part 
of the AGN continuum between them was observed 
through both slits. This produced two independent 
measurements of the 5100—6000 A continuum, since 
the object was centered in each slit separately. 

The spectrograph was rotated to include, together 
with the nucleus of NGC 4151, a field star ("star 1" 
of Penston, Penston, & Sandage 1971, P.A. = 156.3°) 
that served as a local standard. This technique, of 
using a local comparison star, is described in detail 
by Maoz et al. (1990) and Maoz et al. (1994), and 
produces high relative spectrophotometric accuracy. 

For each of the B-side and R-side setups, we tried 
to obtain two or more consecutive exposures per night, 
with typical integration times of 15 min per expo- 
sure. The data were reduced using standard IRAF 
routines with an extraction window of ~ 13". The 
NGC 4151/star ratios of consecutive exposures were 
compared to test for systematic errors. The ratios al- 
most always reproduced to 0.2%— 1.2% at all wave- 
lengths. The few nights where the ratios between ex- 
posures differed by more than 3.5% were discarded. 
(Only three nights were with ratios between 1.2% and 
3.5%.) The NGC 4151 and the star spectra from 
the consecutive exposures were co-added to improve 
the signal-to-noise ratio, and the resultant NGC 4151 
spectrum was divided by an absorption-line-free, inter- 
polated and smoothed version of the co-added stellar 
spectrum. This removes atmospheric absorption from 
the Seyfert spectrum and provides the relative flux cal- 
ibration. This procedure resulted in one spectrum per 
setup per night. Absolute photometric calibration was 
achieved by multiplying the NGC 4151/star ratio by 
a smoothed flux-calibrated spectrum of the compari- 
son star, obtained on photometric nights using known 
spectrophotometric standards. 

The second data set was obtained by the Ohio State 
University (OSU) group with the Perkins 1.8m tele- 
scope at the Lowell Observatory. NGC 4151 was ob- 
served for 10 nights between 1993 December 2 and 11, 
UT, with the Boiler & Chivens spectrograph through a 
5 arcsec slit at P. A. =90°. Many exposures were taken 
each night with integration times of 2 — 3 min each, 
with a total integration time of several hours (in order 
to study the variations on timescales of minutes). The 
spectra covered the wavelength range ~ 4480—5660 A. 



Table 1 
Flux scale factors 





Point- Source 


Extended Source 


Data Set 


Scale Factor 


Correction 


OSU 


0.9558 ± 0.0027 


-0.758 ± 0.018 


Lick 


0.9349 ± 0.0043 


-1.103 ± 0.023 


Wise 


1 






'^In units of 10 ergs s ^ cm ^ A 



The data were reduced using standard IRAF routines. 

The individual spectra were co-added to produce three 
spectra for each night, which, after checking for agree- 
ment of 1% or less, were added into one spectrum per 
night. 

The third data set that will be presented here was 
obtained with the Kast double spectrograph (Miller & 
Stone 1993) on the Shane 3 m reflector at Lick Obser- 
vatory, using a 4 arcsec slit which was aligned along the 
parallactic angle. This set consists of 3 epochs. Seven 
other data sets were obtained by other AGN watch 
members covering anywhere between 1 and 8 nights. 
These will be discussed in the following section. 

2.2. Calibration 

A project of this scope requires the intercalibration 
of the various data sets into a single consistent set. The 
method used for this intercalibration is based on the 
[O III] AA4959, 5007 narrow emission lines, by requiring 
all spectra to have the same flux in these lines. There 
are various ways of achieving this objective. The one 
adopted here is based on Peterson et al. (1994). 

The Wise data set, being the largest one, was cho- 
sen as the reference data set. The [O ill] AA4959, 5007 
emission lines fluxes were measured between the ob- 
served wavelengths of 4952—5051 A. by summing the 
measured flux above a straight line passing between 
the average flux in the 6 A interval on each side of this 
wavelength range. This kind of measurement was de- 
termined by Baribaud & Alloin (1990) to be optimal 
for the [O III] emission lines. 

A second scaling method used is described by van 
Groningen & Wanders (1992). The method flnds the 
optimum scaling factor, wavelength shift, and convolu- 
tion factor of one spectrum with respect to a reference 
spectrum, by slowly varying these parameters until the 
residuals of one or more non-variable narrow lines in 
the difference between the two spectra are minimized. 
This resiilted in spectra and light curves which con- 
flrmed but did not improve (in terms of the internal 
scatter within one data set) the results of the former 
method. 

The mean flux of the two [Oiii] emission lines 
through the Wise apertures was found to be T[o iii] = 
1.575 X 10"-^^ ergs cm~^, with a scatter of 1.1%. 
(This is consistent with the Antonucci & Cohen (1983) 



result of (1.19 ± 0.06) x 10"^^ ergs s"^ cm'^ for 
the [Oiii] A5007 flux.) We consider this to be the 
accuracy of the comparison-star method. As a test, 
all B-side spectra were scaled so that their measured 
J^[o III] agree with the above flux. No obvious differ- 
ences appear in the Wise light curves before and after 
the [O III] scaling, demonstrating the reliability of the 
comparison-star method. An attempt to scale all R- 
side spectra to the same [Sii] AA6716, 6731 emission 
line fluxes gave similar results. A comparison of the 
overlapping continuum bands between the R-side and 
the B-side spectra shows an agreement of about 1% 
for almost all epochs. 

All other data sets were scaled to match their [O iii] 
line fluxes to the Wise ^[oiii] by multiplying each 

J~[OIII] 



spectrum by a factor 



where F{[OIII])ob 



F{[OIII])obs 

is the measured [O ill] line flux of the spectrum. Light 
curves were then measured for the H/3 line and the 
5100—5150 A continuum. These two light curves, 
which were obtained for each data set, were then in- 
tercalibrated with the Wise data set. 

The intercalibration method (Peterson et al. 1994) 
is as follows. First, a point-source correction factor (fi 
was deflned by the equation 



F{HI3) = <pF{HpUs , 



(1) 



where F{Hj3)obs is the H/3 line flux measured for each 
spectrum after scaling its [O ill] line flux to agree with 
the Wise T[o iii], and F(7J/3) is the Wise H/3 flux 
from an observation that is close in time (see below). 
The ip factor accounts for the fact that different aper- 
tures resulted in different amounts of light loss for 
the given point-spread function (which describes the 
surface-brightness distribution of both the broad lines 
and the AGN c;ontiuuum source) and the partially ex- 
tended narrow-line region. 

An additive correction was applied to allow for the 

different amounts of starlight admitted by different 
apertures. This correction, G, was deflned by the 
equation 



Fa(5100) = ^Fa(5100)o6s-G , 



(2) 



where Fx{hlQQ)obs is the continuum measured in the 
observed wavelength range 5100—5150 A after the 



spectrum was scaled to have the Wise [O ill] hne flux, 
and F\(5100) is the contemporaneous Wise resuh. 

A problem unique to this monitoring project is that 
NGC 4151 could be observed by ground-based tele- 
scopes for only about two hours at the end of each 
night. (Due to ROSAT constraints, the campaign 
on this object, which has a right ascension of 12 
hours, was executed in December.) Thus, no pairs of 
nearly simultaneous observations from different obser- 
vatories (to determine the values of ip and G) could be 
found. This limits our ability to intercalibrate the var- 
ious data sets and is different from previous monitor- 
ing campaigns (e.g., NGC 5548) where objects could 
be observed all night, and the time resolution was a 
few days. 

Given the above difficulty, we have taken the alter- 
native approach of interpolation. The reference Wise 
light curves of H/3 and F\(5100) were interpolated lin- 
early, and for every spectrum in the other data sets the 
values of F{Hf3) and i^A(5100) at their precise date 
and time were found, and used for defining and G 
(from eqs. [| and ||) . The ip and G of each data set were 
then determined by averaging the values that were ob- 
tained for each spectrum of a data set. Finally, each 
spectrum of the given data set was multiplied by the 
average ip and subtracted by G, producing the final 
calibrated set. 

The intercalibration constants, ip and G, for each 
data set are given in Tabled. The ip for the OSU data 
set is in good agreement with its value given by Pe- 
terson et al. (1995), who calculated it from the [Olll] 
A5007 surface-brightness distribution of NGC4151 
through different apertures. 

The intercalibration method can only be used for 
data sets with more than two epochs. Based on this 
criterion, several data sets had to be discarded. Also, 
all data sets with an internal accuracy worse than 3% 
(i.e., with spectra from the same night differing from 
each other, or from the Wise set, by more than 3%) 
were discarded. Given the small variation of the source 
during the period of intensive monitoring (1993 De- 
cember 2—11), these data sets do not add useful infor- 
mation and increase the overall noise. 

In view of all these limitations, only the two large 
data sets from Wise and OSU, as well as the smaller 
Lick data set, will be discussed in the following sec- 
tions. 

2.3. Line and Continuum Light Curves 

Fig. I^a and Table ^ give mean continuum light 
curves for four spectral bands: (a) 4560— 4640 A (here- 
after 4600 A), (b) 5100-5150 A (average of the Wise 
B-side and R-side, and several points from the OSU 
and Lick sets, hereafter 5125 A), (c) 6170-6230 A 
(hereafter 6200 A), and (d) 6900-6950 A (hereafter 
6925 A). All wavelengths are in the observed refer- 
ence frame. All four continuum light curves show the 
same feature of small "bumps" superposed on a grad- 
ual rise of the flux throughout the whole monitoring 
period. Examination of these light curves shows that 
the optical continuum changed (from a minimum at 
JD=2449315.6 to a maximum at JD=2449357.6) by 



about 35% at 4600 A and 17% at 6925 A (^^ - 1). 
This can also be seen in Fig. || (bottom panel): the 
r.m.s. of the Wise light curves is ~ 10% at 4600 A 
whereas at 6925 A it is ~ 5%. 

The H/3 flux was measured from each frame of the 
Wise B-side, OSU, and Lick sets, between the observed 
wavelengths 4780-4950 A, by summing the measured 
flux above a straight line passing between the average 
flux in the intervals 4560-4640 A and 5100-5150 A. 
The Ha emission line was measured from the Wise R- 
side and Lick spectra between 6415— 6710 A, and its 
underlying continuum was measured between the av- 
erage flux of two wavelength ranges 6170-6230 A and 
6890-6970 A. Fig. ^ and Table | show the light 
curve, and Fig. |l|c and Table || show the Ha light curve. 
Both Ha and H/3 show a gradual rise in flux through 
the monitoring period, which amounts to '^13% in Ha 
and to - 30% in H/3. While in the H^ light curve 
there are several features which resemble the contin- 
uum light curve (Fig. |^a), such features are barely vis- 
ible in the Ha light curve. 

3. Analysis 

3.1. Cross Correlations 

A main purpose of this and other monitoring cam- 
paigns is to measure the dimensions of the gas dis- 
tribution in the BLR of Seyfert galaxies. This is 
done by cross correlating the line and continuum light 
curves and determining the time lag between them. 
Another goal is to study the variability properties of 
the sources. One of the methods we have applied to 
our data is the cross-correlation algorithm suggested 
by Gaskell & Peterson (1987). In this method, the 
cross-correlation function (CCF) is calculated twice for 
the two observed light curves a{ti) and h{ti): once by 
pairing the observed a[ti) with the interpolated value 
b{ti — r), and once by pairing the observed b{ti) with 
the interpolated value a{ti—T). The final CCF is taken 
to be the average of these two. No extrapolation was 
used to avoid introducing artificial data, and the two 
last observed points of each light curve were omitted 
because of the large separation between them and the 
rest of the light curve. Linear and spline interpolation 
gave similar results. 

A major disadvantage of such interpolation meth- 
ods is the lack of rigorous error estimates for the CGF 
and the deduced lag. One way to estimate the error is 
by assuming a certain BLR geometry and using sim- 
ulations to find the significance of the measured lag 
(Maoz & Netzer 1989). Another way, suggested by 
Gaskell & Peterson (1987), is an analytic estimate for 
the uncertainty on the CCF peak position. This is only 
a rough estimate which relies on specific assumptions, 
such as uniform sampling of the data. We use it for 
lack of a better method. 

An alternative way that avoids interpolation, to find 
the CCF and the time lag with error estimates, is to 
use the Discrete Correlation Function (DCF; Edelson 
& Krolik 1988). The error estimate of this algorithm 
has been questioned by several authors (e.g.. White & 
Peterson 1994; Paper IV). An improved algorithm 



Table 2 
Light curves 



JD 4600 A 5125 A 6200 A 6925 A Ha H/3 
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* OSU data. ^ Lick data. 

^Observation Julian Date — 2449000, rounded to a tenth of a day. 
^Continuum flux in units of 10"^** erg cm~^ s~^ A~^. 
^Emission line flux in units of 10~^^ erg cm~^ s~^. 




Fig. 1. — (a) Four optical continuum light curves for NGC 4151. Wavelength bands are given in upper left corner 
(in A). The ordinate fx is given in absolute units. No shifts are applied to the top three light curves, and the 6925 A 
light curve is shifted by -5xl0~^^ erg cm~^ s~^ A~^. (b) H/3 light curve, (c) Ha light curve. 
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Fig. 2. — Top panel: Mean spectrum of NGC 4151 from the Wise data set. The dotted lines are the mean + 
and — the r.m.s. about the mean. Residual telluric features are visible at 6280 and 6860 A. The continuum and 
line measuring bins are marked. Bottom panel: The smoothed ratio of the r.m.s. and the mean spectrum. The 
variation of the light curves's amplitude and its wavelength dependence can be seen from the larger variance at blue 
wavelengths. 
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Fig. 3. — CCFs of two continuum wavelength bands: 
4600 A and 5125 A, with linear interpolation (solid 
line), spline interpolation (dashed line), and ZDCF (er- 
ror bars with filled circles). Note the good agreement 
between the three methods. 

was recently suggested by Alexander (1996). This new 
approach applies Fisher's z transformation to the cor- 
relation coefRcient, and bins the DCF by equal pop- 
ulation bins rather than equal time bins. It results 
in a more robust and statistically reliable method, the 
Z-transformed Discrete Correlation Function (ZDCF). 
The ZDCF peak position and its errors are estimated 
by a maximum likelihood method that takes into ac- 
count the uncertainty in the ZDCF points. 

Discrete binning implicitly assumes that the spacing 
between the data points is uncorrelated with their ob- 
serving times. The NGC 4151 data treated here pose a 
special problem in this respect. The galaxy's sky posi- 
tion during this project allowed it to be observed only 
at the end of the night. As a result, most of the JDs of 
the Wise set are at 0.6 of a day, whereas most JDs of 
the OSU set are at 0.0 day (Table |). Since the OSU 
data only extend over one eighth of the monitoring pe- 
riod, time lags of n -I- 0.6 days (where n is an integer), 
which cross correlate OSU with Wise points, strongly 
depend on the behavior of the light curve during the 
OSU period, whereas time lags of n days, which cross 
correlate Wise with Wise points and OSU with OSU 
points, reflect the overall behavior of the light curve. 
This special sampling pattern resulted in spurious fluc- 
tuations between consecutive ZDCF points. Such fluc- 
tuations disappeared only when the ZDCF bin size was 
enlarged, or when the OSU data were omitted from 
the light curves. Since the first option significantly de- 
creases the number of ZDCF points, we use only the 
Wise data in the subsequent ZDCF time series analy- 
sis. Below we present results from the ZDCF and the 
interpolated CCF methods (where data points from 
all sets were used). Both methods yield very similar 
results. 

Table ^ shows the properties of the various CCFs 
calculated for the different continuum spectral bands. 
In this table, every wavelength band listed in a col- 
umn is cross correlated with every wavelength band 



listed in a row, such that if the wavelength band in 
the row lags behind the wavelength in the column the 
time lag is positive. For every two bands we calculated 
the CCF using the linear interpolation (column "lin" ) 
and the ZDCF (column "ZDCF" ) methods. For each 
CCF we give the time lag measured from the main 
peak, the error on the time lag (as described above), 
the full width at half maximum of the main peak (mea- 
sured at half way between the peak's maximum and 
minimum), and the correlation coefficient of the peak. 
An example of a CCF of two continuum wavelength 
bands is given in Fig. I The diagonal in Table | gives 
the auto-correlation function (ACF) results for all four 
continuum wavelength bands. The ACFs are plotted 
in Fig. 

Table ^ lists the properties of the CCFs of the 
5125 A continuum light curve with the emission line 
(Ha and II/3) light curves, as well as the CCF of the 
emission lines with each other. This continuum band 
(which was observed in both the B-side and the R-side 
of the Wise set) was chosen since it is the most reliable 
one. 

The Ha light curve (Fig. He) shows a small gradual 
rise of 12% throughout the entire campaign. No 
significant signal was found from cross correlating it 
with itself. In contrast, the H/3 light curve (Fig. |^b) 
shows, on top of the gradual rise of ^ 30%, also some 
features that follow the continuum variations. The 
CCFs of Ha and H/3 with the continuum light curve 
are illustrated in Fig. ^. 

3.2. Continuum Variability 

Our monitoring of NGC 4151 for over two months 
with a temporal resolution of 1 to 4 days resulted in 
light curves with errors of 1%, and can reveal varia- 
tions as low as 3%. The variability timescale, defined 
as the FWHM of the main ACF peak (Table | and 
Fig. Q), is ~ 13 days. The measured continuum vari- 
ations are of 17% to 35%, with amphtude decreasing 
towards longer wavelengths (Fig. |^a). This trend ex- 
tends to shorter wavelengths, as seen by comparing the 
variability at UV wavelengths (Paper I) with the opti- 
cal light curves presented here (Fig. ^). Fig. 0, which 
shows the average and r.m.s. spectrum of NGC 4151 
from the Wise data set, illustrates this variable am- 
plitude (see the larger variance at blue wavelengths in 
the bottom panel). Fig. ^ presents this phenomenon 
in three dimensions: wavelength, JD, and flux. Here 
we have made use of some more continuum wavelength 
bands, and the surface was interpolated over missing 
data points and smoothed. 

To quantify the variability amplitude, and its wave- 
length dependence, we calculated the power spectra of 
the four continuum light curves. Since the first 38 days 
of the monitoring period are evenly and almost regu- 
larly sampled by the Wise data set, it is reasonable to 
apply a discrete Fourier transform to these data. The 
power spectra are presented in Fig. |^. The decreasing 
power of the variability with increasing wavelength is 
clear. A power-law fit, of the form PDS cx to the 
first six points gives an index a of —1.5 ± 0.9 for the 
wavelength bands 4600 A and 5125 A, and 
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See explanation in text. 



Table 4 
Emission line CCFs 
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See explanation in text. 




Fig. 4. — ACFs of the four continuum bands presented in Table ||. Notation is the same as in Fig. |. (a) 4600 A, 
(b) 5125 A, (c) 6200 A, and (d) 6925 A. 
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Fig. 5. — CCFs of the 5125 A continuum band with 
Ha (top panel) and H/? (bottom panel). Notation is 
the same as in Fig. |[ Note the time lag of ^ — 3 
days indicated by the main peak. 



-0.8±0.6 for the wavelength bands 6200 A and 6925 A. 
The small variation amplitude, especially at the longer 
wavelengths, prevents us from determining a signifi- 
cant dependence of variation timescale on wavelength. 

The dependence of the variability amplitude on 
wavelength can be explained by the varying relative 
contribution of starlight from the underlying galaxy. 
We show this by using template spectra of Sab galax- 
ies from Coleman, Wu, & Weedman (1980) and Kin- 
ney et al. (1996). Based on previous observations of 
NGC 4151 (see Maoz et al. 1991) we estimate the 
galaxy contribution to the spectrum through the Wise 
aperture, in the continuum wavelength band of 4600 A, 
to be in the range (2 — 3.5) x 10^^"* erg cm^'^ s~^ A^^. 
These numbers are confirmed by Peterson et al. (1995) 
who find the starlight contribution through the Wise 
aperture to be 2.2 x 10^ erg cm^ -2 A-i. We 
normalize the template spectra to these numbers and 
subtract the appropriate values from the four contin- 
uum light curves. When using the maximum value 
for the galaxy contribution at 4600 A we find that 
the relative change in each light curve during the 14 
days of intensive monitoring is about 20%, i.e., the 
differences in amplitude between the light curves pre- 
sented in Fig. 1^ disappear. The 20% variation ampli- 
tude is also similar to that of the UV continuum at 
2625-2750 A (hereafter 2688 A) but is different from 
that of the 1250-1300 A (hereafter 1275 A) light curve 
where the relative variation is ^35%. 

We illustrate this result in Fig. ^ where we present 
a simulation of the continuum light curves based on 
the 2688 A continuum light curve. From the normal- 
ized galaxy template we separated the total flux in 
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Fig. 6. — Same as in Fig. pL for the period of the 
lUE campaign. Two UV continuum light curves are 
shown for comparison with the optical continuum light 
curves. 



this light curve into a galaxy contribution 2%) and 
AGN contribution (- 98%). The resulting AGN light 
curve was scaled by a factor according to its part in the 
other wavelength ranges and was added to the galaxy 
contribution at each wavelength range. (The upper 
limit of the galaxy contribution at 4600 A, of 3.5 x 
10""'^'' erg cm~2 A~^, gives the solid lines in Fig. O, 
and the lower limit, of 2 x 10^^^ erg cm^^ s^^ A , 
gives the dotted lines.) The observed data at all opti- 
cal wavelengths show good agreement with the galaxy- 
diluted UV light curve. The wavelength dependence 
of the variability amplitude in the range 2700-7200 A 
can thus be explained by the different starlight contri- 
bution at different wavelengths. The 1275 A UV con- 
tinuum light curve's amplitude cannot be explained 
using the starlight contribution alone, as illustrated in 
the top panel of Fig. ^. 

From the CCFs of the different continuum wave- 
length bands with each other, no significant lags are 
detected, i.e., the continuum varies in phase at all opti- 
cal wavelengths, to within 1 day. There is no apparent 
lag between the optical and UV light curves, as dis- 
cussed in 3.3 and in Paper IV. 



No clear conclusions result from comparison of the 
optical light curves to the high-energy light curves (Pa- 
per III), since the high-energy monitoring period was 
brief compared to the optical variation timescale, and 
the number of data points is relatively small. 
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Fig. 8. — Power Density Spectra for the four contin- 
uum wavelength bands (calculated points and power- 
law fits). Note the decreasing power in variability with 
increasing wavelength. 

3.3. Emission Line Variability 

The Ha light curve (Fig. |^c) shows a small gradual 
rise throughout the entire campaign, and its ACF peak 
has a width of the duration of the program. No vari- 
ability timescale can be deduced from this ACF. Cross 
correlation of Ha with the 5125 A continuum gives a 
time lag of - 2 days (Fig. | top panel). The H/3 
ACF indicates a variability timescale similar to that 
of the continuum (Table ||). Its CCF with the 5125 A 
continuum (Fig. ^ bottom panel) suggests a lag, with 
respect to the continuum light curve, of — 3 days. 
The CCF of the Ha and the H/3 light curves shows no 
lag, implying that both emission lines varied in phase. 

We have performed cross correlations of the UV 
continuum light curves (presented in Fig. ^ with the 
Ha and H/3 emission-line light curves. No significant 
time lags were found from these CCFs' peaks. The 
strong variations in the UV continuum light curves are 
of order 1—2 days and such rapid variations were not 
detected in the optical emission-line light curves. This 
suggests a BLR size larger than 1—2 It-days, in which 
light-travel-time effects smear out the rapid variations 
of the ionizing continuum. Since no large timescale 
variations (of order 10 days) took place during the UV 
monitoring period, no time lag between the UV light 
curves and the optical emission lines could be found. 

Detailed analysis and cross correlation of the optical 
light curves with the UV light curves (Paper I), for 
the lUE monitoring period, are shown and discussed 
in Paper IV. Here we discuss only the correlation of 
the UV light curve with the entire optical period that 
include 14 days of monitoring at the Wise Observatory, 
prior to the beginning of the lU E campaign. 

Cross correlation of the entire optical 5125 A contin- 
uum light curve with the UV 1275 A and 2688 A con- 
tinuum and the UV Hues Civ A1549 and Heii A1640 
(Paper I) are shown in Fig. |l^. All three CCFs, of 
2688 A, Civ, and Hen, show the same feature of a 
broad maximum, ranging from to 6 days. We be- 
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Fig. 9. — Light curves from Fig. g superposed on the 
linearly-interpolated 2688 A continuum light curve, af- 
ter it was scaled and diluted with a constant starlight 
flux according to the galaxy contribution to the dif- 
ferent wavelength bands: solid lines correspond to a 
galaxy contribution at 4600 A of 3.5 x 10^^^ erg cm^^ 
s^-'^ A"-'^, and dotted lines correspond to a contribution 
of 2 X 10~^^ erg cm~^ s~^ A~^. Note the agreement 
between the scaled UV light curve and the optically 
observed points (bottom panel), indicating that the 
decrease in variability amplitude at larger wavelengths 
is mainly the result of dilution by the constant galaxy 
light. The 1275 A UV light curve (top panel) is qual- 
itatively different from the scaled 2688 A light curve. 

lieve that this apparent positive delay is an artifact 
due to the finite duration of the monitoring campaign. 
This broad feature is probably a result of the optical- 
continuum behavior prior to the UV monitoring: these 
three UV light curves show a gradual rise followed by 
period of constant flux and the optical-continuum light 
curve shows a gradual rise prior to the beginning of the 
UV monitoring. Those similar shapes result in high 
correlation (~ 0.8) over a period of a few days. 

We have carried out two tests of this hypothesis. 
We have extrapolated the UV 2688 A continuum to 
the period prior to the lUE monitoring based on the 
shape of the optical continuum. This was done by 
finding several 5125 A continuum data points which 
coincide in time with several 2688 A continuum data 
points. The fluxes of those pairs of points were fitted 
with a linear function which was then used to scale the 
5125 A light curve prior to the lUE campaign to the 
2688 A level. This extrapolated UV light curve was 
cross correlated with the optical continuum and with 
the C IV line light curves, and gave a CCF peak at zero 
lag. The second test is cross correlation of the optical 
and UV light curves using only the part of the optical 




Fig. 10 — CCFs of the entire optical continuum 5125 A band with the (a) 1275 A UV continuum, (b) 2688 A UV 
continuum, (c) C IV A1549 hne, and (d) He ii A1640 line. Notations are the same as in Fig. |3[ The spline interpolation 
CCF is not shown. 



light curve which overlaps the lUE campaign. As 
shown and discussed in Paper IV §3.4, the results are 
narrower peaks CCFs consistent with a zero time lag. 
Thus, both tests suggest that the time lags between 
the optical continuum and the UV line and continuum 
light curves are consistent with zero. The last conclu- 
sion explains the narrow peak (0-2 days) CCF of the 
5125 A continuum with the 1275 A continuum shown 
in Fig. |lo|a. The zero time lag is, indeed the lag and 
the difference between this CCF and others shown in 
Fig|l^ are due to the fact that the 1275 A continuum 
has more features in its light curve than the other UV 
light curves. 

3.4. Comparison with Past Results 

Maoz et al. (1991) monitored NGC 4151 in 1988 
for a period of 8 months with an average sampling fre- 
quency of once every 4 days. They found continuum 
variations of ~20%, similar to the 1993 amplitude, 
with a typical timescale of ^^28 days, about twice what 
is found in the present campaign (~13 days). They 
also found that both Ha and H/3 followed the contin- 
uum variations and were able to determine time lags 
of 9±2 days with respect to the optical continuum, for 
the two emission-line light curves. In this work we find 
that the Ha and H/3 variations lag behind the optical 
continuum variations by — 3 days. In both campaigns 
the H/3 response to the continuum variations is larger 
than the Ha response. 

To investigate whether these discrepancies in 
emission-line lags are the result of the different con- 
tinuum variability timescales or of different analysis 
methods, we have applied the auto-correlation meth- 
ods described above to the Maoz et al. (1991) data. 
We confirm that when using the interpolation methods 
(linear and spline) the variation timescale is ~28 days. 
However, using the ZDCF we find a variation timescale 
of ~15 days. This difference can be explained by the 
fact that the 1988 monitoring program had several 
large gaps which, when interpolated, increased the cor- 
relation between the continuum and line light curves 
at larger times. Using the data from the last 96 days of 
the 1988 campaign, which are more regularly sampled 
and have only one large gap, we find that both the 
interpolation and the ZDCF methods yield the same 
variability timescale of ^15 days. NGC 4151 therefore 
had similar optical continuum variability behavior in 
1988 and in 1993. 

We have also performed cross correlation of the 
emission lines and the continuum of the 1988 data. 
All cross-correlation methods yield the same results 
for the time lag, ~ 8 days, consistent with the lag 
found by Maoz et al. (1991). The lag remains the 
same when only the last 96 points of the 1988 data are 
used. We conclude that the differences in Ha and H/3 
lags behind the continuum between 1988 and 1993 are 
not the results of different methods of analysis. 

Finally we have checked by way of simulations 
whether the small lag we measure for the emission lines 
in 1993 could be the result of the particular form of the 
continuum variation in 1993, which consisted mainly 
of a monotonic rise. Assuming a spherical BLR shell 



of inner radius 2 It-days and outer radius 30 It-days 
(found by Maoz et al. to best fit their results), we cal- 
culated the emission-line light curve produced by such 
a geometry when driven by the spline-interpolated 
4600 A continuum light curve from the present cam- 
paign. We then sampled these two light curves at the 
same epochs of our actual observations, and added to 
the sampled points a noise similar to the errors in our 
light curves. We applied to these light curves the var- 
ious cross-correlation techniques discussed above, and 
found the CCFs to have broad peaks at lags of ~ 10 
days. Thus the above thick shell geometry, with an as- 
sumed ionizing continuum similar to the 4600 A light 
curve of 1993, gives a lag which is much larger than 
that which we measure. The small time lags found 
in the 1993 data are probably not an artifact of the 
particular continuum behavior we observed. 

If the above differences in emission-line lag are real, 
they may be related to the different state of NGC 4151 
which, at the time of the 1993 campaign, was in an "ac- 
tive" state, in contrast to its lower flux level during the 
Maoz et al. monitoring campaign of 1988. (The opti- 
cal continuum and line fluxes in the present monitor- 
ing campaign are about a factor of two higher than in 
1988; see also Oknyanskij, Lyutyi, & Chuvaev 1994.) 
Evidence for a changing lag has been shown by Pe- 
terson et al. (1994) in the one other Seyfert galaxy 
that has been intensively monitored for several years, 
NGC 5548. A possible physical explanation for the 
change in the lag is a real change in the BLR gas dis- 
tribution between 1988 and 1993. Such a change is, 
in principle, possible considering the scales and veloc- 
ities present in the nucleus of NGC 4151 which yield 
a dynamical time scale of ~ 3 years (see also Wanders 
1994). An alternative explanation is that during the 
present campaign, the optical continuum did not prop- 
erly represent the behavior of the ionizing continuum. 
In particular, if the variability timescale of the ion- 
izing continuum driving the lines is different between 
the two campaigns, the resulting time lags for the same 
BLR geometry can be very different (Netzer & Maoz 
1990; Netzer 1990). Some evidence for this is seen in 
the differently-shaped light curves of the 1275 A and 
the 2688 A continua (Fig. P) . This possibility is exam- 
ined further in Paper IV. 

4. Summary 

We have presented optical-band results of an inten- 
sive two-month spectrophotometric monitoring cam- 
paign of the Seyfert galaxy NGC 4151, with a typical 
temporal resolution of one day. The main results of 
this campaign are as follows. 

1. The continuum variations are between 17% and 
35%, with decreasing amplitude towards longer wave- 
lengths. The broad Ha line flux varied by 12% and 
the broad H/3 flux by ~ 30%. 

2. The decreasing continuum variability found at 
longer wavelengths can be explained by the varying 
contribution of starlight from the underlying galaxy. 
The exception to this is the far-UV 1275 A continuum 
where the variations must be intrinsically larger. 

3. The various optical continuum bands vary in 



phase, with a lag of < 1 day. The typical continuum 
variability timcscalc is ^ 13 days and is similar at all 
optical wavelength bands. The variability amplitude 
and timescale are similar to those observed in the past 
in this object. 

4. No evidence for a time lag between the optical 
continuum and the UV continuum and emission lines 
was found. This may be partially the result of the 
short duration of the lUE campaign. Paper IV gives 
details on intcrband phase lags derived from the CCFs. 

5. The Ha and H/3 light curves follow roughly the 
continuum variations and lag them by — 3 days, in 
contrast to past results where a time lag of 9±2 days 
was found. This may be related to a different variabil- 
ity timescale of the ionizing continuum, or to a real 
change in the BLR gas distribution in the 5.5 years 
interval between the two campaigns. 
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